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Abstract 



. We derive finite temperature expansions for relativistic fermion systems in the 

00 

■ presence of background magnetic fields, and with nonzero chemical potential. 



We use the imaginary-time formalism for the finite temperature effects, the 
proper-time method for the background field effects, and zeta function regu- 
larization for developing the expansions. We emphasize the essential difference 
between even and odd dimensions, focusing on 2 + 1 and 3-1-1 dimensions. 
We concentrate on the high temperature limit, but we also discuss the T = 



^ ■ limit with nonzero chemical potential. 

I. INTRODUCTION 

The study of fermion systems in the presence of external electromagnetic fields has ap- 
plications in diverse areas of physics, including astrophysics, solid state, condensed matter, 
plasma and particle physics. Indeed, astrophysical considerations led to the first systematic 
study of relativistic noninteracting fermion systems The existence of very high inten- 
sity magnetic fields in gravitationally collapsed objects motivated further investigations of 
the energy-momentum tensor and equation of state for a degenerate electron gas in a strong 
uniform magnetic field . More recently, this problem has been addressed in the framework 
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of finite temperature quantum field tlieory , using the fact that many of the thermody- 
namic properties may be derived from the corresponding finite temperature effective action. 

In this paper we present a pedagogical discussion of the finite temperature (T) and finite 
chemical potential (/i) effective action for fermions in an external static magnetic field. Our 
aim is to review and unify the imaginary-time approach developed initially for finite 

temperature free systems , with the Fock-Schwinger proper-time approach ||9|-[TTll developed 
for systems interacting with external electromagnetic fields but at zero temperature. For- 
mally, these two approaches fit together beautifully for the case of static external magnetic 
fields; in fact, the finite T and finite /x effects are completely separate from the computation 
required to compute the effects of the external static magnetic field. Thus, knowledge of the 
zero T and zero fi proper-time effective action is completely sufficient to write down a corre- 
sponding expression at nonzero temperature and nonzero chemical potential. However, this 
expression for the effective action is formal and the separation between T and fi effects and 
those of the external field may become blurred when one tries to make various approximate 
asymptotic expansions, such as for high or low temperature. Such expansions are necessarily 
complicated due to the proliferation of energy scales: thermal energy kT, fermion mass m, 
chemical potential fi, cyclotron energy B/m, and also possible momentum scales associated 
with spatial variations of the external field. We concentrate mainly on the high tempera- 
ture limit, in which T (we use units in which Boltzmann's constant k = 1) is the dominant 
energy scale. However, we also present a simple approach to the other extreme: the T = 
limit. This complements Refs. which have focussed on the low T limit, and which 

have primarily used the real-time formalism for discussing the finite temperature effects. 

Throughout our analysis, we treat both 2 + 1 dimensional and 3 + 1 dimensional theories. 
This is motivated by the known profound differences between 2 + 1 and 3 + 1 dimensions 
(and in general between odd and even dimensional space-times) for free fermion systems 



at finite T and /i [0,|T^, and for zero T fermions in external fields ||TT|. We find that a 
consistent treatment of both cases requires careful use of zeta function regularization, which 
has been used previously for free systems [|l2|,|l3l (and for systems with constant external 



Ao [0). For the magnetic backgrounds we also find that the 2 + 1 and 3 + 1 cases involve 



very different expansions at high temperature. The high temperature behavior of QED2_,_i 
is also of interest for studying questions of spontaneous symmetry breaking |T5HI9| . 



In Section II we review the structure of the zero temperature effective action with nonzero 
chemical potential for 2 + 1 and 3 + 1 dimensional fermions in external static magnetic 
fields. Finite temperature is introduced in Section III using the imaginary-time formalism. 
In Section IV this is combined with the proper-time formalism to provide a general formal 
expression for the finite T and n effective action. In Section V we apply this to the high 
temperature limit of the free fermionic theories, and in Section VI to the high temperature 
limit of fermions in a static magnetic field. The zero temperature limit is examined in 
Section VII, and we conclude in Section VIII with some comments regarding possible further 
extensions of this approach. 

II. EFFECTIVE ACTION FOR MAGNETIC SYSTEMS 

The basic object of interest in this paper is the effective action 

z^eff = LogDet (z^-m-/i7°) (1) 

where p= Dy'^'^ = [dy + leA^) 7^, and we choose Minkowski gamma matrices 7*^ satisfying 

{7^ 7^^} = -2g^'^ = 2 diag(l, -1, -1, ... , -1) (2) 

The term /i7° in the Dirac operator in (|l]) reflects the presence of a chemical potential /i, 
corresponding to a term —^ip'^ip in the Lagrangian. 

In 2 + 1 dimensions, the irreducible gamma matrices may be chosen to be 

7° = a3 ^^=ta^ 7^ = ^^^ (3) 

where the a* are the 2x2 Pauli matrices. Note that an alternative choice, 7° = — o"^, 
7^ = — zcr^, 7^ = —ia"^, corresponds to changing the sign of the mass, m —>■ — m, in the Dirac 
operator appearing in the effective action (|Tp. The system with effective action (|I|) is not 



parity invariant since a fermion mass term breaks parity in 2+1 dimensions |jT5|-|T^]. However, 



a parity invariant model may be constructed by considering two species of fermions, one of 



(4) 



mass m and the other of mass —m (see Footnote 11 in Ref. [|T5|). This may be achieved by 
choosing a reducible set of 4 x 4 gamma matrices 

'Y 
-Y 

in which case the Dirac operator is block diagonal: 

ip — m — fij^ 

—ip — m + ;U7° 

Thus, the effective action for this parity invariant system may be written as 



(5) 



iSes = Log Det2+i [iD^T^ — m — /iF^ 



LogDet2+i - (iDo - fiy + + [D ■ ^ 



(6) 



In 3 + 1 dimensions parity symmetry is not an issue, but the effective action (|l]), which 
involves the first-order Dirac operator, may still be written in the same form as (^, which 
involves a second order operator. We use the fact that there exists an additional gamma 
matrix 7^ satisfying {7'',7^} = and (7^)^ = 1. Then 

Det3+i {tp-m- /i7°) = Detg+i (7' {tp-m- /i7°) 7^) = Detg+i {-tp- m + /i7°) (7) 

Therefore, the effective action (|l|) may be expressed as 



^•S'eflf = -LogDeta+i (ip- m - /i7°] (-i^ 

= ilogDeta+i (- {tDo - fif + + ■ 7) 



m + fi'y 
2 



2 
2 

Note that the spatial operator (^D ■ 7^ which appears in (|^) and in 
operator. In 2 + 1 dimensions it reduces to 

D_D+ 
D+D^ 

where the magnetic field is B = F12 = diA2 — d2Ai, and 



(8) 



is a positive 



(5 



•7 



(9) 



D± = Di± tD2 

In 3 + 1 dimensions the operator (^D ■ 7^ reduces to 
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(10) 



(z3.7) 



(11) 



If we choose the external magnetic field to be directed along the direction and to be 
independent of x^, then (with a suitable choice of gamma matrices) this may be simplified 
further to 

\ 



-dl- 



( D_D+ 

D+D_ 
D_D+ 

V D+D_ ) 



(12) 



where D± are as defined for the 2 + 1 dimensional system in ([ToD . 

Thus, in each case, the spectrum of the spatial operator [d ■ 7^ is determined by the 
spectrum of the 2-dimensional Schrodinger-like operators D±D^. The operator [d ■ 7^ is 
effectively diagonal and we may write 



(13) 



For static magnetic backgrounds = and we can replace {iDq — jj) in (|^) or (|^) by — /i, 
where uo is an energy eigenvalue. Therefore 



iSpff = I -— Tr Lo, 



27T 



(14) 



Here, the trace operation Tr is understood to mean a summation over the eigenvalues 
of both the positive operators — D^D^ and — D^D+ in 2 + 1 dimensions, resp. 



— 9? — D,D_ and — fi? — D_D+ in 3 + 1 dimensions. 



III. FINITE TEMPERATURE FORMULATION 

It is clear from ( [T^ ) that the effects of the external static magnetic field are contained 
solely within and are clearly separated from the chemical potential yU and the energy 
trace over uj. We can therefore pass to a finite temperature formulation just as in the free 
case 0-0, by replacing the energy integration with a discrete summation: 

duj_ i_ ^ 

27ri / 1\ 



Here 13 = 1/T, where T is the temperature and Boltzmann's constant k has been absorbed 
into T. The transition to finite temperature would not be so straightforward if there were 
external electric fields, but here we consider only external static magnetic fields. Also note 
that in the zero temperature expression (|1^) it looks as though the dependence on the 
chemical potential fi may be formally eliminated through a naive shift of the integration 
variable uj. However, such a shift would violate the boundary conditions used to compute 
the trace, and a proper treatment at zero temperature leads to the appearance of non- 
analytic behavior in fi, corresponding to sharp cut-offs in the energy spectrum ||7|j20|]. At 
finite temperature, these sharp cut-offs are smoothed out, and the dependence on /i is 
correspondingly smooth, as we shall see below. 

The effective action ([1^) may now be expressed as 

'5'eff = ^ J2 TrLog -{ujn+ {\S\ - fl)^{uJn- {\S\ + fi)^ 

P n=—oo 
-| oo 

-TrLog n [((1^1 - + + l)/pf) ((|^| + ^^f + (vr(2n + l)/pf) 

n=0 

-Tr Log 



1, 



cosh — (l^^l — /i) cosh — {\E\ + fj.) 
2 2 



(16) 



where in the last step we have used the infinite product representation of the cosh function 



cosh(x) 



n 1 

n=0 \ 



(17) 



and we have dropped an infinite contribution that is independent of S"^. 



It is a simple matter to re- write (16) as (dropping again an irrelevant constant) 



S, 



eff 



Tr 



1^1 + ^Log (l + e-^(l^l-^)) + ^Log (l + e-'^(l^l+^)) 



(18) 



This expression for the effective action generalizes an analogous expression (with /i = 0) 
derived in [^] for the free case. The only effect of the chemical potential is to shift the 
'energy' eigenvalue \£\ by =i=/i, which corresponds to a shift in the threshold energies for 
particles and antiparticles. The only effect of the external static magnetic field is to modify 
the spectrum of the 'energy' eigenvalue \£\ from the free spectrum to a spectrum involving 
dependence on the external B field. Therefore, at least in principle, we can now compute 
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the effective action for any external static magnetic field for which we know the spectrum 
of the operators D±D^ appearing in (^) and (|I2|). 

When fi = and we take the zero temperature limit {f3 — oo) then the effective action 



in ( ]T8| ) reduces to 

^eff U=0;/3^oo =Tr(|£|) (19) 

which is the familiar T = 0, fi = effective action in a static magnetic background [ p!0[] . 

When fi and j3 ^ oo we must distinguish between the cases fi < m and fi > m. 
When fi < m, all low temperature thermal excitations are exponentially suppressed because 
|£^| ± /i > (since \S\ > m). Therefore, in this case the effective action reduces just as 
in (|19|). However, if /x > m, then in the infinite (3 limit the first logarithmic term in (|T8[) 
contributes for the portion of the spectrum for which \S\ — fi < 0. Thus, we have 

'S'cffU^0;/3-.oo = Tr \S\ + {fi-\S\)9{fi-\S\) 

= tt fi e{fi - \s\) + \s\ e{\s\ - fi) (20) 

where 6 is the step function. This is the standard expression for the effective action at 
zero temperature and with nonzero chemical potential [^0- The first equality in (pOD 
emphasizes the correction from the zero T and zero fi expression (p!9D, while the second 
emphasizes the physical content of the effective action with zero T and nonzero fi as fi times 
the number of occupied particle states plus the trace of the energy eigenvalues above the 
threshold fi ||20|. In Section VII we examine this T = limit in detail for fermions in a static 
magnetic background. 

The form of the effective action ([T8| ) is reminiscent of the grand partition function in 
non-relativistic statistical mechanics. Indeed, the non-relativistic limit corresponds to the 
situation in which the rest mass energy m is the dominant contribution to \S\, 



\£\ = \lm'^ +(D-if = m + ^(D-^Y + ... (21) 



In this limit, with fi> m and m — >■ cxd, the first logarithmic term in ([18|) dominates over the 
second {i.e antiparticles are suppressed), and we are left with 



^eff ^ -^Tr Log (^1 + e^'^N«e-^(^-^) j (22) 

where we have identified /xnr = yU — as the non-relativistic chemical potentiaL The expres- 
sion (^2]) is times the logarithm of the grand partition function for the corresponding 
non-relativistic fermion system. 



IV. PROPER TIME FORMULATION 



The proper time formulation provides an efficient method for computing the effective 
action at zero temperature PHTII, and furthermore has the virtue that the generalization 
to finite temperature and nonzero chemical potential naturally separates out the infiuence 
of a static background magnetic field. Using an integral representation of the logarithm to 
define the logarithm of an operator, we may express the finite temperature version of the 
effective action (|T^ as 



S, 



cS 



oo 

- TrLog 



n=— oo 
oo 



1 ^ f^ds 



exp s -s 



(23) 



where we have chosen to refer all energy scales to m in order to have dimensionless operators 
in the exponent. For massless theories one must choose a different reference energy scale, as 
discussed in Section VI. 

Expression (p3D may be re-cast in terms of the second elliptic theta function 92{u\t), 

Airis 



cS 



1 f°° ds f 27r/is 



Tr 



exp 



(24) 



where 21,22 



^2(M|r) = 2 ^ ei-r(n+l/2)2 ^^3^ + l)u) 



(25) 



ra=0 



A Poisson summation formula converts the second theta function into a fourth theta function 
according to the identity 



u 



-1/2 



(26) 



This converts the expression 

5'cfr = 



for the effective action into a form involving 64,: 



m 

2^ s 



°° ds 1 1/3 fi 



Tr 



exp 



1^1 



(27) 



where the fourth theta function is 



dA{u\T) = l + 2^(-l)V"™' cos(2nu) 



(28) 



n=l 



This expression (^) for the effective action is particulary useful as it shows clearly the 
separation between the effects of the static background magnetic field, which appear solely 
in the trace factor, and the effects of finite temperature and nonzero chemical potential, 
which appear solely in the 64 '^J^ ^ factor. Indeed, when /i = and /9 — > 00, the 64 
factor reduces to 1 and we are left with the standard proper time expression for the zero 
temperature effective action in a static system. This corresponds to keeping just the term 



"1" in the expansion ( pS]) of the 64 function, so the remaining summation over n in (pHf ) 
represents the nonzero temperature correction. The utility of elliptic theta functions in the 
computation of finite temperature effective actions has been noted previously in for 
fermions (without chemical potential) and in |^ for bosonic systems. 

All information about the static magnetic background is neatly encapsulated in the 
proper time propagator 



Tr 



1^1 



exp 



Tr 



exp 



(D . 7)^ 



(29) 



which is computed independent of any reference to temperature or chemical potential. Thus, 
if one computes the zero temperature effective action using the proper time method it is 
completely straightforward to then write down an expression for the effective action at finite 
temperature and at nonzero chemical potential simply by inserting the 64 factor as in (pTj). 

However, while the expression ( pT] ) illustrates the separate roles of finite /3, fi and the 
external static field, it is not so straightforward to use it to obtain useful numerical estimates. 
This is because of the wildly oscillatory behavior of the 6'4-function in ( p7\ ) for large values 
of the proper time parameter s. This oscillatory behavior is also sensitive to the magnitude 
of the dimensionless parameter which appears in the first argument of the 64 function. 
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These difficulties are further comphcated by the prohferation of energy scales - for zero T 
and /i the only scales are the fermion mass m, the characteristic strength scale B (with 
dimensions of m^) of the external magnetic field, and possibly also characteristic length 
scales associated with spatial variations in the magnetic field. The generalization to nonzero 
temperature and chemical potential introduces two further energy scales: (3 and fi. In this 
paper, we concentrate mainly on high temperature expansions (in which T = l/j3 is the 
dominant energy scale), although we also discuss the zero temperature limit (with nonzero 
chemical potential) in Section VII. 

To conclude this brief review of the finite temperature formalism for fermionic systems, 
we show how the general expression ( pT]) relates to the previous expression (0), for which 
our statistical mechanics intuition is most direct. Using (p8|) and rescaling the proper time 
variable s in (0), s/m^ —>■ jS'^s/A, we obtain 

1 r°° ds ' 



S, 



off 



/3v^ Jo s3/2 

1 r°° ds 



9, 







g-/32£-2s/4" 


\ 2 


TTS I 





1 + 2 £(-l)"e-"'/^ cosh(n/5/i) 



n=l 



Tr 



-/32f2s/4 



(30) 



The integrations over s may now be performed, with those in the summation term requiring 



the identity (see H Eqs. 3.471.9 and 8.469.3) 



°° ds 



n Jo s 



3/2 



exp 



n 
s 



-I3n\£\ 



n> 



n 



(31) 



Thus 



S-cff = Tr|^| + - ^ e-^"l^l cosh(n/?/i) 



which is just expression ([TSl). 



(32) 



V. FREE THEORIES 

Before discussing the temperature expansions for fermions in the presence of background 
magnetic fields, we first describe the high temperature expansions for free theories. This 
is partly to establish some notation and to introduce some number theoretic functions (the 

10 



zeta and eta functions), but also in order to point out some important subtleties that arise 
even in the free theories and which have implications for the magnetic systems. 



In the free fermionic theory in {d + l)-dimensional spacetime, 



m 



+ p ^ , where p i 



IS 



a (i-dimensional momentum vector. Thus, the proper-time propagator ( p9D is simply 



Tr 



exp -s 



c f m Y 



(33) 



2 \2^/^ J 

where c is the spinor dimension (the d = 2 and d = 3 cases discussed previously correspond 
thus to c = 4). Therefore, the effective action (pTf ) is 



S, 



eff 



m \ r°° ds 

V^j Jo 7¥ 



1 + 2 ^(-1)" cosh(/3/in) exp 



n=l 



4s 



(34) 



This expression already illustrates an important difference between the 3 + 1 and 2 + 1 
dimensional cases, for which (c? + 3)/2 is an integer or half-odd- integer, respectively. For 
example, in the zero temperature and zero chemical potential case, the square parentheses 
in (0) reduce to a single term "1" , and so the 3 + 1 free case naively leads to a divergent 
factor r(— 2) which must be regulated consistently. This may be achieved, for example, by 
cutting off the lower limit of the s integration, or by shifting the dimension c? ^ 3 + e 
and extracting a finite piece. In contrast, in 2 + 1 dimensions the corresponding factor is 
r(— 3/2), which is finite. These issues are well understood for the zero temperature theories 
PHTTI, but below we illustrate analogous differences between d odd and d even for the (3 and 
11 dependent contribution to the effective action. 

The s integrations in the n summation in (0) may be performed using the following 
integral representation of the modified Bessel function (see [^, Eq. 7.12.23) 



TT 



Re z > 0, arg^; < — 
> I 6 I 2 



(35) 



This leads to|^ 



^As mentioned previously, the 1 in the square bracket in Eq. (34) corresponds to the zero tem- 
perature (and zero chemical potential) proper time representation of the effective action, which we 
denote by S^^^'^~^ . 
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Sr-ff — S, 



r=0;At=0 



cflf 



2c 



m 



d+l oo 

I El 



-1)" cosh( J3 fin) Kd+i (Pmn) 



' f3mn^ 



-{d+l)/2 



(36) 



For free theories the only energy scales are m, fi and T = so a high temperature 

expansion corresponds to (3m <^ 1 and jSfi <^ 1. Therefore, expression (^) may be converted 
into a high temperature expansion by using the ascending series expansions of the modified 
Bessel function Kiy{z). The difference between 3 + 1 and 2 + 1 dimensions is reflected 
in the fact that the series expansions of K,y{z) are very different for i/ an integer oi v a, 
half-o dd-int eger . 

For V = integer = N, the expansion of Ki^{z) for small z begins with a z~'^ term and 
also involves a logarithmic piece: 



Kr,(z) 



1 f 

2 V 2 / „_n j\ \ 4 / 



j=0 



-1 



N oo 
i=0 



j!(iV + j)! 



log I - ^^{j + 1) - ^V^(Ar + J + 1) 



(37) 



where ^p{x) = ^logr(x) is the digamma function ||2T|. On the other hand, for u 
half-odd-integer = + |, the corresponding small z expansion is 



K 



N+ 



2 hj^-jy- 2^- 



1 oo 

-y 

2 ^ 



r(iv-j + i) fz\-^-y^ r(-iv-j 



2 



2\ J 



(38) 



which begins with but has no logarithmic piece. 

These Bessel function expansions (^7|) and (^Sj) lead to the following explicit high tem- 
perature expansions of the free effective action (^6]). In 2 + 1 dimensions, or more generally 
in (i + 1 dimensions with d even. 



od even qT=0;/^=0 

•^eff — ^eff 



m 



d+l 



5:(-l)"cosh(/?/in) 

.2V^/ n=l 



oo 



3=0 



r(^ 



j) / /?mn 



'^'^'^ , r(-(^)-j) 



2^2^2\ J 



13 m n 



(39) 



j! V 2 ; J! 
Notice that the cosh(/5/i?T,) term can be expanded in powers of n^, in which case the sum- 
mation over n has the form of the eta function |^ (an alternating series analogue of the 
Riemann zeta function) 
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OO -| 

^W = E(-ir"'z7 (40) 

n=l " 

For Rer > 1 the eta function //(r) is related to the Riemann zeta function, ({r) = J2'^=i ^/n^, 
by 

r/(r) = (1 - 2i-X(r) (41) 

For negative integer values of r, this function takes especially simple values. For even 
negative integers 

7^(0) = ^, v{-^k)=0 k = 1,2,3,... (42) 

while for odd-integer negative arguments 

r^(l - 2k) = {2'' - 1)^ k = 1,2,3,... (43) 

where is the {2k)*^ Bernoulli number. Also note that ri{l) = log 2. 

Since ri{—2j) vanishes for j = 1,2,..., only the j = term survives in the second 
summation term inside the square brackets in (^9]), and this surviving term then cancels 
against the zero T, zero /z effective action. In contrast, the first summation term inside the 
square brackets in (^) involves odd powers of n, and so contributes for all j. Finally, we 
obtain the following double expansion 



-<d even 
'efE 



-id+1) ^ (/5/i)2fc (-lyri^ -j)7]{d + i -2k- 2j) ((3-^ 



which agrees with the corresponding expansion in Eq. (4.24) of Ref. [|T^], although the sum 
is organized in a different manner. Notice that the leading term at high temperature is 
proportional to T'^^^, with the remaining terms in the expansion involving higher powers of 
j3m or (3^. For zero chemical potential the expansion simplifies to a single sum 



-(d+i) ~ (-l)Jr(!i±i - i)r]{d + 1 - 2j) Il3 



Sir" = c Y. ''" "^ (45) 

and in the massless case 

^cff \m=0 = C [fWT^) r( — ) 2^ j^^^ (46) 
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In 3 + 1 dimensions, or more generally for d+1 dimensions with d odd, the Bessel function 
expansion ( |37D leads to the following high temperature expansion of the free effective action: 



S^,""' = S:,='-^^=' -2c[-^] E (-1)" cosh(/?/.n) 



m 



2^ 



X 



n=l 

d-l 



2I 2 j j! I, 4 j 

r i)(d+i)/2 (/3'm^nV4)^ Aog ^ - -W;r7 + 1) - -th(n + 



(47) 



The first term inside the square brackets produces T'^^^ times a polynomial Pd+iiPfi, (3m) of 
order d+1 in and /3m, because the expansion involves only even powers of ra, for which the 
eta function vanishes if these are positive even powers, but gives a finite contribution while 
the power is negative or zero. Thus the summation over j and the expansion of cosh(/9/in) 
are each truncated. The effective action takes then the form 



+ 



Tayl(/?m,/3/i) [ (48) 



r((d + 3)/2) 

where the Taylor series Tayl(/3m, /3/i) is determined using the following further summation 
identity (for details see [1T2| ): 

.og(2.„) ^ { ^^l^ _ i«)Wc(2. + 1). («) 



VI. CONSTANT MAGNETIC BACKGROUND 



As discussed in Section in order to obtain a closed expression of the form (^) for 
the finite temperature and finite chemical potential effective action for fermions in a static 
background magnetic field one simply needs the proper time propagator (pQ]). For the case 
of a uniform background field this is a well-known computation whose exact solubility is 
due to the fact that the spectra of the operators —D±D^ correspond to one-dimensional 
harmonic oscillators, with 'frequency' 25, each energy level having degeneracy given by the 
(infinite) background magnetic flux ||2^ . 
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- = 2Bn 

-D^D+ = 2Bn 

Thus, the proper-time propagator is 



n = 1, 2, 3, . 



n = 0,1,2,3, 



(50) 



Tr 



exp -s 



(51) 



3 + 1 dim. 

f coth (If) , 2 + 1 dim. 

The extra (2^) "TS factor in 3+1 dimensions is due to the integration over k^, the momentum 
corresponding to the free motion in the z direction and 3/271 is the usual Landau level 



degeneracy factor |26 



The effective action (l27| ) is therefore 

\ d-l 



eflf 



/ m ^ 



B 



ds 



271 Jo s('^+l)/2 



, , Bs\ im? 
coth 1 — - ^ 
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4:71S 



(52) 



where d = 2 or d = 3. Note that we have explicitly subtracted the zero B contribution 
by subtracting the leading small s divergence m?/{Bs) from the function coth(i?s/m^). 
This ensures that the effective action vanishes with vanishing 5, and simply corresponds to 
computing a difference of the effective action with and without the B field. We could also 
subtract the B^ contribution as this may be absorbed into the original action via a charge 
renormalization. However, we choose to defer this additional subtraction for two reasons. 
First, we wish to emphasize that the B^ contribution differs in an interesting manner for 
d = 2 and d = 3. Second, it is conventional [P-[TT[| to make this subtraction only for the 
zero T and zero /i portion of the effective action, since in the 3 + 1 dimensional case this 
contribution involves a divergence whose regularization involves a charge renormalization. 

To develop high temperature expansions from (|52|) , along the lines discussed for the free 
case in the previous Section, would require expanding the 64 factor in the integrand and then 
performing the proper time integrals. However, these integrals cannot be done in closed form 



(in contrast to the free case where they led to modified Bessel functions - see (|35|) ) because 



of the coth.{B s / m"^) factor. We therefore rescale the proper time parameter s: s/m? 
Then in the high T limit we expand the coth(/5^i?s) factor: 



coth B'Bs 



(3^Bs 



00 r)2k r2 
\ - ^ 02k 



ph. 



(53) 
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where B2k are the Bernoulh numbers This leads to the following double expansion 

fc=i y^'^')- n=i ^ [fJmny 2 

Assuming that the temperature T dominates all other energy scales, we can further expand 

using the series representations (|37D and (^) of the modified Bessel functions K^{z). Note 

that, just as in the free case, there is a significant difference between d = 2 and = 3, as the 

former case involves Bessel functions with half-odd-integer index, while the latter involves 

Bessel functions with integer index. 

For d = 2, 

e.=2 _ 4 ^ B,,{^r ^ {-^y r (1 - J - 2A:) ^ ip^y^r^jS - 2j -4k- 21) 

Notice that the zero T, zero fi effective action is canceled by a term arising from the expansion 
of the Bessel functions. Considered as an expansion in the magnetic field, the leading term 
is of the form 

Stf = (3B'fiPf^,Pm) + --- (56) 

where f{f3^,(3m) is a dimensionless function whose leading term is a (negative) constant. 
The expression (^5|) has a well-defined massless limit: 

4 .fi B,Ul3'B/2rr (I - 2k) ^ (i3^fr,(3 - 4k - 21) 

And for zero chemical potential, 

_ 4 ^ B,,(l3'B/2r fi (-gf )' r (I - J - 2k) „(3 - 2, ^ 4fc) 

For d = 3 the situation is somewhat different owing to the appearance of a divergent 
term in the zero T, zero /i effective action and of logarithmic terms in both the zero T, zero 
piece and in the finite temperature piece of the effective action. From (|5^) , using the 
Bessel function expansion (^7|), we find 



T{Ak + 2j + 21-3) 



X (l - 2^'=+2,+2/^3 ^^4^ + 2j + 2/ - 3) 



r(2A; + j 
16 



(59) 



except the term 
fc^l j^O Z-0 



Notice that the T = 0, /i = effective action cancels against a term 



nr ^ 



(2A:)! 



(60) 



arising from the Bessel function expansions, leaving a term quadratic in the B field with an 
(infinite) coefficient 

1 



C 



dss e 



1 „~m 



(61) 
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which may be absorbed in a renormalization of the electric charge P]. 

The remaining expansion (|59|) has no quartic or quadratic term in the temperature but 
has a logarithmic dependence on the temperature. This agrees with Ref. but disagrees 



with Ref. 

To conclude this discussion of the high temperature limit for fermions in a constant 
magnetic background, we observe that in the 2 + 1 dimensional case it is possible to give an 
alternative derivation of the magnetic background effective action beginning instead with 
the expression (p^B]). For ease of presentation, consider the case /i = 0. Then 



^ 2 /3 2^(22^^-1)^ 

= -TrLog cosh-|£:| = (^^^(^p)! ^ 



|2p 



(62) 



The trace over the even powers of the energy operator can be expressed in terms of gener- 
alized zeta functions. In 2 + 1 dimensions 



Tr 



2p 



B 
2^ 



2Y,{m'^ + 2Bn)P + m 



2p 



n=l 



m 



2Ti M ^' 2B 



B 
2^' 



-m 



2p 



(63) 



where the generalized zeta function is defined as |^ 



c(^,?) = E 



(64) 



^0 + 

When the first argument of the generalized zeta function is a negative integer, it reduces to a 
polynomial (in the second argument) which is in fact proportional to a Bernoulli polynomial: 

Bp+i{x) 



The Bernoulli polynomials are given by 



p+l 



(65) 
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Bp{x) 



p / p 

E 



(66) 



r=0 \r 

with expansion coefficients involving the Bernoulh numbers Br- 
Thus, the high temperature expansion ( p2| ) may be written as 

2 \ , 1 / 2\ P' 

p + 1 1 m"^\ 



Qd=2 



eff \fJ.=0 



1 g (22P - 



2p 



p+1 



i3 



p+i 



m 
2B 



+ 



2B 



(67) 



2vr/3%tlP(P+l)(2p)! 

Remarkably, this expansion (|67D is simply another way of expressing the double expansion 
in (|58|) , up to a 5 independent term. 



VII. ZERO TEMPERATURE LIMIT 

Having considered the high temperature behavior of the effective action, we now turn to 
the other extreme: T = 0. With zero chemical potential this is just the usual field theoretic 



effective action, which may be computed using the proper-time formalism p|-p!l|. With 
nonzero chemical potential, the zero T limit is more subtle due to nonanalytic contributions 
from the sharp cut-off of the energy spectrum which is no longer smoothed out by thermal 
fluctuations. The low T effective action with nonzero chemical potential and background 
magnetic field in 3 + 1 dimensions has been studied recently ^j^, primarily using the real- 
time formalism. Here we discuss a somewhat different (and more direct) approach to the 
low temperature limit, based on the expression (^). According to (|20|) , the correction AS* 
to the zero T and zero fi effective action when the chemical potential is nonzero may be 
written as 



AS\t=o;ii^o = Tr 



- \S\)9{fi - \£\) 



(68) 



Since \S\ > m it is clear that this correction to the usual zero temperature effective action 
vanishes unless ^ > m. For the remainder of this Section we assume that /i is indeed greater 
than m. 

The free case is well known (see for example [^) and we turn directly to the case of a 
magnetic background. We consider the 2 + 1 and 3 + 1 dimensional cases separately. In 2 + 1 
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dimensions, to compute the correction term (^) we simply fill the Landau levels, according 
to the Landau level spectrum (|50D, up to just below the level of the chemical potential 

2B 



AS- 



■d=2 



T=0;^^o = ^{2 J2 (/X - Vm2 + 2Bn) + (/i 



m] 



Bm J /i 



TT m 



/i^ — m? 
2B 

jj, — m 



+ 



I2B J 1 



2'~ ' 2B 

2' 



2B 



2m 



'25 7 1 m 
^^\2'^^2B 



(69) 



The last two terms are in fact canceled by the zero chemical potential part of the effective 



action. We note the appearance of a nonanalytic dependence on the integer part 



2B 



of 



2 2 

the dimensionless ratio ^ ~^ . In the nonrelativistic limit, this leads to oscillatory behavior 



of the effective action and related thermodynamics quantities (such as the magnetization) 
as the strength B of the external magnetic field is varied. 

Specifically, in the nonrelativistic, weak field limit, ^ (//^ — m^) ^ B, the 2 + 1 
dimensional correction term is 



ASf=^ 



5 I ^ ^ J / B 



(70) 



where /xnr = fi — m is the nonrelativistic chemical potential. This may be re-expressed as 



^'S'NR^|T=0;At^O 



B^ 
Tim 





"m/iNR" 




[ B \ 



B^ 
vrm 



B 



^2 oo sin 



m/iNR 



B 



/m/iNR _ 1 
\ B 2 



B 

m/iNR\ 



5/i 



NR 



27r 

-B/^NR 



(71a) 



47r^m 



fc=i 



Txm \ B J 



2ti 



2n 



(71b) 



where we have used the zeta function identity (see |21[] 1.10 Eq. (6)) 



C(s,i;) = 2(27r)"-^r(l-s)^F-^sin(27rfci; + 7rs/2), 0<i;<l,Res<0 (72) 

k=l 

Thus the effective action includes a term which as a function of 1/5 is periodic with period 



1/ (m/iNR), the characteristic period of the de Haas - van Alphen oscillations |23]. It is 
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interesting to notice that (in contrast to 3+1 dimensions) there is an exphcit algebraic 
expression ( |71a] ) for the effective action that shows these oscillations. 
In 3 + 1 dimensions, the relativistic correction term ( p8| ) is 

2B 



AS"'-'! 



27r 



(/" - yP^ + m2 + 2Bn) 



5 



(27r)2 

2B 



jj, + jj? ~ m? 



m 



'/i + \/ — — 2Bn 



(73) 



n=i \ \ yAn^T2Bn ; ; ) 

This expression agrees with a corresponding one in |Q, derived from a zero temperature 
limit in the real-time formalism. Expanding the logarithmic term in the weak field limit (in 
which B <^ ij?) we obtain 

Bm? I yU l^'^ — m? //i + V/i^ — m?' 



AS' 



d=3\ 



\T=o■,^J.^o 



(2 



TT r m V ^ 



- Log 



7? 



7r2 f^^ (2/ + l)(2/ + 3) 2' 25 



m 



/i^ — 



25 



3 /i^ — 
2' 25 



(74) 



This expression involves a term depending on the nonanalytic fractional part of ^ . Again 
oscillations arise when we represent the first zeta function with the help of the zeta function 
identity (^ 



3 /j^-m^ 
2' 2B 



23 



sin 27rA; , „ 
2r(/ + 5/2)^ ^ ^ 

k=l 



2B 



Til 37r 

2 4 



(75) 



(27rA;)'+5/2 

These oscillatory components reveal themselves in a nonrelativistic limit ^ (/i^ — m^) 3> 
B where only the first term / = survives, being responsible for the usual nonrelativistic de 
Haas - van Alphen effect p7[] : 



d=3 loscill 



NR IT=0;m^O - 4^^4 



55/2 00 sin (27rA;^™ - f ) 



E 

k=l 



J^5/2 



(76) 
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VIII. CONCLUSIONS 



In this paper we have shown how to compute finite temperature expansions of the effective 
action for fermions in the presence of a static background magnetic field. The computational 
techniques of the imaginary time formalism (for the finite temperature effects) and the proper 
time method (for the background field effects) combine elegantly to produce a simple formal 
expression (^) for the finite T, and nonzero chemical potential, effective action. Explicit 
high (and low) temperature expansions may then be obtained given the explicit functional 
form of the proper time propagator (pUj). Generically, high temperature expansions involve 
logarithmic terms in d + 1 dimensions with d odd, with these terms being absent if d is 
even. These logarithmic terms may be traced to corresponding terms appearing in the series 
expansions of the modified Bessel function Kj\f{z). For a uniform magnetic field background, 
the zero temperature limit (with nonzero chemical potential) may also be evaluated in a 
straightforward manner in terms of generalized zeta functions. These zero temperature 
effective actions involve oscillatory pieces which reflect the Landau level structure of the 
underlying single particle spectrum. 

The formalism developed here is motivated in part by the fact that the proper time 
method provides a powerful technique for analyzing non-constant background fields. For 
example, one could generalize from a uniform static background magnetic field to a static 
magnetic field with spatial dependence. Such a generalization has been considered recently 
for the zero temperature (and zero chemical potential) effective action in 2 + 1 dimensions, 
with one particular nonuniform background field leading to an exact closed-form expression 
for the effective action [^. It would be interesting to incorporate the finite temperature 
and finite chemical potential effects into such a model. Since expression (^Tj) completely 
separates the finite temperature (and nonzero chemical potential) effects from those of the 
static background, it is very tempting to make high (and low) temperature expansions at 
this level, without having an explicit expression for the proper time propagator. However, 
even in the free case this leads to subtleties which can be missed if one is not careful. For 
example, even in the free case, one might think to make a high temperature expansion by 
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expanding the exponential term inside the proper time integral in the effective action (p4D. 
This amounts to deriving the series expansions (|37| ) and (^8]) for the modified Bessel function 
K^{z) directly from the integral representation p5| ) by expanding in powers of z"^ inside the s 
integral. But this type of expansion does not produce the necessary logarithmic terms when 
u is an integer - instead, to obtain the correct series expansion one must use a different 
integral representation based on a contour integral [^. For fermions in the presence of 
some general background gauge field, such naive expansions, made inside the proper time 
integral, must also be treated with care, particularly in 3 + 1 dimensions. 
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